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ABSTRACT 

Long-term orbital evolution of multi-planet systems under tidal dissipation often converges to a 
stationary state, known as the tidal fixed point. The fixed point is characterized by a lack of oscillations 
in the eccentricities and apsidal alignment among the orbits. Quantitatively, the nature of the fixed 
point is dictated by mutual interactions among the planets as well as non-Keplerian effects. We 
show that if a roughly coplanar system hosts a hot, sub-Saturn mass planet, and is tidally relaxed, 
separation of planet-planet interactions and non-Keplerian effects in the equations of motion leads to 
a direct determination of the true masses of the planets. Consequently, a "snap-shot" observational 
determination of the orbital state resolves the sin(J) degeneracy, and opens up a direct avenue towards 
identification of the true lowest-mass exo-planets detected. We present an approximate, as well as a 
general, mathematical framework for computation of the line of sight inclination of secular systems, 
and apply our models illustratively to the 61 Vir system. We conclude by discussing the observability 
of planetary systems to which our method is applicable and we set our analysis into a broader context 
by presenting a current summary of the various possibilities for determining the physical properties 
of planets from observations of their orbital states. 

Subject headings: planets and satellites: general — celestial mechanics — methods: analytical 



1. INTRODUCTION 

Since the seminal discovery of the first giant planet or- 
biting a main sequence star (Mayor & Queloz 1995), us- 
ing the radial velocity (RV) method, over 400 additional 
extra-solar planets have been confirmed. The greatest 
disadvantage of the RV method lies in the uncertainty of 
the true masses of the discovered planets, as the inclina- 
tion of the orbits to the line of sight, /, are unknown. 
In resonant systems, such as GL876, monitoring of the 
resonant argument and the precession rates may lead to 
determination of the true masses (e.g. Rivera et al 2005). 
In the vast majority of cases, however, the sin(I) degen- 
eracy, remains a continued source of frustration. 

Still, RV surveys persist in yielding fruitful results, and 
the continued detection of exo-planets has brought forth 
many surprises. Perhaps one of the biggest surprises has 
been the discovery of extremely close-in bodies whose 
mass-range spans the entire planetary spectrum. These 
objects have since become a subject of fascination in the 
community and more importantly, have provided a new 
test-bed for various theoretical efforts. 

Extra-solar multi-planet systems that host "hot" plan- 
ets differ drastically from our own solar system in many 
ways, including orbital dynamics. In our solar system, 
gravitational interactions among the planets are suffi- 
cient to, at least approximately, explain orbital evolu- 
tion. In many extra-solar planetary systems however, 
similarly to the case of the Galilean satellites, dissipa- 
tion of orbital energy due to tides plays an unavoidably 
important role. The long-term effect of this additional 
interaction provides an opportunity to infer important 
additional properties of the system that cannot be ob- 
served directly. 

Qualitatively speaking, in a system of two or more 

kbatygin@gps.caltech.edu 



planets that are not in a mean-motion resonance and 
are roughly coplanar, tides drive the orbits towards a 
stationary state i.e. a "fixed point" . A fixed point is 
characterized by continued apsidal alignment and a well- 
determined eccentricity ratio that is nearly constant in 
time (Wu & Goldreich 2002, Mardling 2007). The fac- 
tors that determine the actual quantitative nature of the 
state are not limited to gravitational planet-planet inter- 
actions. Indeed, general relativistic and tidal corrections, 
among other things, play a crucial role. It is through 
these "non-Keplerian" interactions that additional infor- 
mation can be learned, as they are governed by parame- 
ters other than just planetary masses,. 

Upon discovery of the first multiple planetary system 
with a transiting "hot Jupiter", Hat-P-13 (Bakos et al 

2009) , it was pointed out that the system likely resides at 
a fixed point (Batygin, Bodenheimer & Laughlin 2009). 
Furthermore, it was shown that as the mass and radius 
of the inner planet are known, consideration of the plane- 
tary quadru-pole gravitational field, and its contribution 
in determination of the fixed point leads to a direct mea- 
surement of the planetary interior structure. In other 
words, a precise "snap-shot" of the orbits gives the plan- 
etary Love number, fc 2 , which is a measure of the interior 
density distribution, with high accuracy. 

The last decade of observations has revealed that gen- 
erally, hot Jupiters tend not to be accompanied by read- 
ily detectable companion planets (Ragozzine & Holman 

2010) . Smaller planets, such as hot Neptunes and hot 
Super-Earths, however, tend to occur in multiple-planet 
systems (Lo Curto et al 2010), hinting at different mi- 
gration histories (Terquem & Papaloizou 2007). 

Here, we consider the latter class of systems, with an 
eye toward inferring conventionally unobservable plan- 
etary properties that influence the details of the fixed- 
point configuration. In particular, we show that if a 



2 



Batygin & Laughlin 



non-transitng (RV) system hosts a small (R < -^Nep) 
hot planet, it is possible to derive the true masses of the 
planets i.e. resolve the sin(J) degeneracy from a detailed 
determination of the system's orbital state. The plan 
of the paper is as follows: in section 2, we outline our 
mathematical model. In section 3, we apply the theory 
to the 61 Vir (Vogt et al 2010) system. In section 4, we 
discuss the possibility of determination of the radius and 
interior structure of massive RV planets. We conclude 
and discuss our results in section 5. 

2. DYNAMICAL EVOLUTION OF A PLANETARY SYSTEM 
WITH A CLOSE-IN PLANET 

As already mentioned above, there are important dif- 
ferences between the dynamics of systems with and with- 
out close-in planets. In conservative (Hamiltonian) sys- 
tems, of which our solar system provides an excellent 
approximation (Laskar 2008), the motion of the planets 
is subject to Liouville's theorem. Accordingly, strictly 
Hamiltonian flow can have no attractors in phase-space 
(Morbidelli 2002). Conversely, in dissipative systems, the 
phase space volume explored by the system continuously 
contracts, and truly steady-state solutions are possible. 
In other words, tides are needed for the system to arrive 
to a stationary state. 

Thc path that the system will take to the fixed point is 
non-unique and depends on the initial conditions. Con- 
sequently, the initial transient period will also depend on 
the initial state. However, the fixed-point itself is unique 
for a chosen set of system parameters, and the system 
has no memory of its own evolution once it arrives to 
the fixed point. Thus, any quantity that is inferred from 
the fixed point is independent of the system's formation 
history. We now describe a mathematical model for the 
system's evolution to a stationary state and its orbital 
characteristics. 

2.1. Secular Interactions With non-Keplerian Effects 

Whenever planets are far away from low-order mean 
motion commensurabilities and the orbits are not chang- 
ing significantly on the orbital time-scale (i.e. planets are 
not scattering), a secular approximation to the dynam- 
ics can be made. The secular approximation refers to an 
averaging procedure, where the gravitational potential 
between planets is averaged over the mean longitudes, 
thereby reducing the degrees of freedom inherent to the 
problem. 

Since the pioneering work of Laplace (1772) and La- 
grange (1776), a number of perturbation theories based 
on various approximation of the disturbing potentials 
have been developed and applied in both solar system 
and exoplanetary contexts (Le Verrier 1856, Brouwer & 
van Wocrkom 1950, Laskar 1986, Laskar 2008, Eggleton 
& Kiseleva-Eggleton 2001, Mardling & Lin 2002, Lee & 
Peale 2003, Michtchenko & Malhotra 2004, Migaszewski 
& Gozdziewski 2009, Lovis et al 2010, etc). Still, it is 
perhaps easiest to illustrate the ideas presented here in 
the context of a modified Laplace-Lagrange (LL) secular 
theory. 

The classical secular disturbing function (planet-planet 
potential), of N secondaris that interact solely by Newto- 
nian gravity, expanded to first order in masses and second 



order in eccentricities reads (Murray & Dcrmott 1999) 
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where e is eccentricty, w is the longitude of perihelion, a 
is semi-major axes and n is mean motion. The constant 
coefficients A take the form 
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where aj k = aj/a k if (aj < a k ); a k /aj if (a k < aj), 
otj k = a jk if (aj < a k ); 1 if (a k < aj), b^) 2 (a jk ) & 

b 3 / 2 (aj k ) are Laplace coefficients of first and second kind 

respectively, and m — rh/ sin(J) are the true masses of 
the planets i.e. rh are the measured minimum masses 
and I is the inclination of the system from line of sight. 

Upon application of thc linear form of Lagrange's 
planetary equations in terms of polar coordinates (h = 
ecos(w),k = esin(tt7)), a linear system of ODE's 
emerges, where the A matrix encapsulates the dynamics 
of the system: 

fc=i y " fe=i 

We can express the system of equations more compactly 
by switching to complex Poincare variables z = e e ra = 
h + ik. Simple chain rule yields 

dz N 

~£ =^J A jkZj k (5) 

This eigensystem can be solved in the standard way, sim- 
ilar to the problem of N coupled pendulums, and the 
solution reads: 

N 

z j {t)=Y J ^e i{aM) (6) 
k 

where g's are the eigenfrequencies and /J's are the eigen- 
vectors of the A matrix. The relative amplitudes of the 
eigenvectors and the corresponding phases, 6, are deter- 
mined by initial conditions. 

The above formulation does not take into account the 
additional orbital precession induced by general relativ- 
ity (GR), stellar and planetary spin, and the tidal bulges 
of the star and thc planet. The classical LL solution of- 
ten gives poor quantitative approximations to the orbital 
evolution of extra-solar planets (Veras & Armitage 2007), 
where the additional precession can dominate (Ragozzine 
& Wolf 2009). 

The contributions to apsidal precession from the 
above-mentioned effects can be written as follows (Sterne 
1939): 

fdvj\ 3GM*n 
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where c is the speed of light, k 2 is the Love number (twice 
the apsidal motion constant), R is physical radius and 
f2 is the spin frequency. In equations (8) and (9) the 
first terms correspond to the planet and the latter terms 
correspond to the star. Neglecting higher-order effects, 
the total additional apsidal precession, evaluated for each 
planet can be organized into a square diagonal matrix 
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and added to the A matrix in equation (5). This matrix 
is not to be confused with the mutual inclination interac- 
tion matrix (see Murray & Dermott 1999), for which the 
standard notation is the same. When evaluating the ad- 
ditional precessions (the B matrix) in the context of LL 
theory, it is customary to expand equations (7) - (9) to 
first order in e, such that the dependence on e disappears, 
and equation (5) remains linear in eccentricity, thus re- 
taining its analytical solution. The augmentation of the 
diagonal matrix coefficients will modify the eigensystem 
quantitatively. However, the qualitative essence of the 
solution remains unchanged: the solution (equation 6) is 
still a sum of sinusoids with constant amplitudes. 

So far, we have retained all additional precession terms 
for the sake of completeness. Before proceeding further, 
let us examine the relative importance of the terms that 
depend on the physical properties of the bodies with re- 
spect to GR, which is a purely geometrical effect and 
only depends on stellar mass and the orbital parameters. 
Consider the following dimension-less numbers 
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A Jupiter-like planet at a characteristic close-in orbit 
(P ~ 3 days) has A p spm ~ 0.05 A* spin - 3 x IO" 5 , 

^ 5 x I0~ 4 . Inflated hot Jupiters 
3> I, due to the R p dependence 
of the tidal term. Thus, precession rates of many hot 
Jupiters are completely dominated by the planetary tidal 
term, distantly followed by GR. As mentioned above al- 
ready, this effect has has been used to infer the interior 
structures of transiting hot Jupiters, both in isolation 
(Ragozzine & Wolf 2009) and in the presence of a per- 
turbing companion (Batygin, Bodenheimer & Laughlin 
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2009). Conversely, for a Neptune-like planet on a 3-day 
0.005 A* pin ~ 3 x IO- 5 , A p udal ~ 0.1 and 

3 x 10~ 6 . The numbers continue to decline for 
super-Earths and terrestrial planets. This implies that 
in practice, the apsidal advance, resulting from rotation 
of both the planet and the star, as well as that result- 
ing from the stellar tidal bulge, can often be neglected. 
Indeed, the situation is bimodal: for large planets, tidal 
precession dominates, where as for small planets, GR 
dominates the extra apsidal advance. 

Let us now add dissipative tides to the system. Gen- 
erally, tidal heating conserves the total angular momen- 
tum, but not energy. This leads to decay of the planetary 
eccentricity, as well as decay (or growth, depending on 
stellar spin) of the planet's semi-major axis (Goldreich 
1963). The evolution of the semi-major axis happens 
over a much longer time-scale than that of the eccentric- 
ity, so in our simplified model, we adopt the standard 
practice of parameterizing tides with a constant decay of 
the eccentricity, dz/dt = z/t c , where t c is the circular- 
ization timescale (Goldreich & Soter 1966): 
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Here, P is the orbital period and Q is a tidal quality 
factor. In a similar fashion as above, each planet can 
be subjected to tidal damping of eccentricity by con- 
structing a square diagonal matrix with the elements 
Cjj = 1/t^\ Note that because tidal dissipation only 
affects semi-major axes, eccentricities and rotation rates 
directly, an identical procedure cannot be carried out for 
the mutual inclination eigenmode solution (see Mardling 
2010 for an in-depth discussion). The equation of mo- 
tion that accounts for the additional precession and tidal 
damping of the eccentricity takes the form: 
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At this point, we have changed the solution qualitatively. 
The introduction of eccentricity damping has added a 
complex component to the eigenfrequencies. Conse- 
quently, in the secular solution (6), real exponential de- 
cay factors appear in front of the oscillatory solution. 
The eigenvectors are now damped. Furthermore, the 
imaginary components of the eigenfrequencies need not 
be equal, and generally will not be, except for a narrow 
set of system parameters. This implies that the decay 
timescale of one of the modes, 



(Imfe])~ 



(15) 



can be considerably longer than all others, and the sys- 
tem will eventually evolve to a state that is characterized 
by a single eigenmode. Note that the eigenmode de- 
cay timescale can greatly exceed the tidal circularization 
timescale, prolonging the lifetime of the dissipated plan- 
ets' eccentricities. Upon inspection of equation (6), it is 
clear that once the system is characterized by a single 
eigenmode, the rates of orbital precession are identical 
for all planets in the system. From Lagrange's planetary 
equations, this automatically implies that the apsidal an- 
gles between the orbits must be equal to Aw — or 
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Fig. 1. — A damped, modified Laplace-Lagrange secular solution 
of a 2-planet system with mi = 10 _5 MQ,m2 = 10 -2 Mq, ai = 
0.03AU, a 2 = 0.3AU (a = 0.1) and e 2 = 0.1. Three solutions are 
presented corresponding to the initial conditions e\ = 0.1 (Blue), 
e\ = 0.03 (Red), and e\ = (Green), with randomly chosen lon- 
gitudes of perihelia. The black line shows the eccentricity of the 
outer planet. The apsidal angles initially circulate, but switch to 
libration at t fa Q X 6 X 10 J years. The system reaches a fixed 
point as the anti-aligned (gi) mode decays away completely at 
t Ri Q X 1.3 X 10 6 years. Note that the system looses memory of its 
initial conditions as it approaches the fixed point. 




Fig. 2. — Fixed point eccentricity ratio as a function of sys- 
tem inclination for a 2-planet system with m\ = 10 -5 MQ,m2 = 
lO _2 M ,ai = 0.03AU and a 2 = 0.3AU (a = 0.1) (see Fig. 
1). The black curve, labeled LL was computed directly from the 
Laplace-Lagrange eigenvector solution. The blue dashed curve is 
the approximation to the LL solution, given by equation (19), cor- 
responding to e 2 = 0.1. The curves with eccentricity labels demon- 
strate the dependence of the eccentricity ratios on the stationary 
eccentricity of the outer secondary, as dictated by the secular per- 
turbation theory, developed by Mardling (2007). Recall that m 
refers to the RV minimum mass. 



Aw = it. In other words, all orbits are either aligned or 
anti-aligned, depending on which particular eigenmode 
has survived. Additionally, in this case, the ratios of the 
eccentricities are also well-defined by the eigenvector of 
the surviving mode. When the system has reached a state 
where its dynamics are characterized by a single mode, 
it has reached a "fixed point." Addition of higher-order 
terms to the disturbing function will modify the the ec- 
centricity ratios implied by the fixed point, but will not 
cause Am to be anything other than or n. 

2.2. Determination of sin(I) 

We now have all the necessary ingredients to determine 
the system inclination. Usually, the dissipation time- 
scale greatly exceeds the secular time-scale (C <C A,B), 
so a system at a fixed point is characterized by a sin- 
gle eigenvector of the [A + B] matrisQ From the defini- 
tions of the coefficients of A (equations 2 & 3) , it is clear 
that they are linearly proportional to sin(/). In fact, we 
can replace the true masses, m, by the minimum masses 
m in equations (2 & 3) and write A — A/ sin(I). B 
is however independent of the system inclination, given 
that GR is the only contributing factor. Recall that this 
is the case for Neptune-sized and smaller planets, for 
which A spin <C 1 and A ti( j Q ; <C 1. As a result, the eigen- 



vectors of the 
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matrix, which physically 



A/sin(I) 

correspond to the eccentricity ratios of the planets, de- 
pend explicitly on sin(/). Namely, every value of the 
system inclination corresponds to an eccentricity ratio of 
the planets. Consequently, a precise observational deter- 
mination of the eccentricity ratios yields the true masses 
of the system. Let us turn to an illustrative example 
below. 

2.3. Beyond Linear Order in e: 

the Case of Well- Separated Orbits 



1 The physical effect of including C in the solution is to offset 
the apses by a small factor proportional to Q~ 1 . 



Consider the case of two well-separated (a -C 1) sec- 
ondaries, where the inner planet is on a close-in orbit. In 
such a scenario, we only need to consider the additional 
apsidal precession of the inner planet. Since a -C 1, it is 
sensible to expand Laplace coefficients in equations (2-3) 
into hypergeometric series and retain only the first terms: 
b^ha) « 3a, bfUa) « (15/4)a 2 . With a little algebra, 
it is easy to show that to leading order in a and 77, the 
eigenfrequencies take on a simple form: 
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where F = Bu/An and 77 = Cn/An. The two eigen- 
frequencies physically correspond to modes dominated 
by the inner (<?i) and outer ((72) apsidal precessions. 
Note that the imaginary components of the modes have 
explicitly different dependences on a. The multiplier 
in equation (16) is just An, expanded to first order 
in a. So neglecting sin(J) for the moment, it is clear 
that Im [gi] = l/r c . This is consistent with the ob- 
servation of Mardling (2007) that ~ 3t c are needed 
for the system to attain a stationary state. The situ- 
ation is wildly different however for the second mode, 
aslm[p 2 ] = (25/16)(m 1 /m 2 )(a 5 / 2 /r c )(l + T)- 2 . Con- 
sequently, equation (15) implies that t 
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unless mi 3> TTI2 and the overall lifetime of the inner 
eccentricity is also greatly enhanced. 

The corresponding eigenvectors, also to leading order 
in a, but neglecting the higher-order correction from 77 
read: 
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Note that the eigenvector of the first mode is negative. 
By Euler's identity, the negative sign introduces an ad- 
ditional iir in the exponent of the solution (6) for one of 
the planets. Physically, this corresponds to apsidal anti- 
alignment. Thus it is apparent from equations (18 & 19) 
that the first and the second eigenmodes correspond to 
anti-aligned and aligned orbits respectively. 

As an illustration, consider a pair of planets with 
masses mi = 10~ 5 Af©/sin(I), m 2 = 10~ 2 M Q /sin(J), 
and semi-major axes a\ = 0.03AU, a 2 — 0.3 orbiting 
a M ± = IMq star. The (/ = 0) damped, modified 
Laplace-Lagrange secular solution of this system is pre- 
sented in figure (1), where the planets were started with 
ei = e 2 = 0.1 and randomly chosen longitudes of per- 
ihelia. The planetary Love number was chosen to be 
k 2p = 0.3. Let us examine the evolution in some detail. 
After an initial transient period of ~ 3t c , the system 
reaches a fixed point. Thereafter, the free eccentricity 
decays on the timescale of r = Im \g?\ s» 10 5 t c . As al- 
ready stated, the addition of a perturbing planet has 
prolonged the lifetime of the dissipated planet's eccen- 
tricity immensely. As a result, it must be pointed out 
that the detection of an eccentric close-in planet alone 
does not imply that the planet itself is weakly dissipa- 
tive. Rather, self-consistent calculations are required to 
place any constraints on Q. 

The above analysis implies that planets on well- 
separated orbits in a tidally relaxed system will be ap- 
sidally aligned rather than anti-aligned, with the fixed- 
point eccentricity ratio, ei/e2, given by the correspond- 
ing eigenvector. Figure (2) shows the solution for the 
eccentricity ratio as a function of system inclination, /. 
The solid line, labeled LL, represents the directly cal- 
culated eigenvector and the dashed line represents the 
approximate solution, given by equation (19). 

As can be inferred from figure (2), and equation (19), 
the fixed point eccentricity of the inner planet is much 
smaller than that of the outer planet. This is troublesome 
in the context of LL theory, where the outer eccentricity 
is already assumed to be small, because a precise obser- 
vational determination of the eccentricity ratio becomes 
difficult. Consequently, we need to lift the constraint 
on the outer secondaries' eccentricity, so that the inner 
one at least becomes observably large. This can be ac- 
complished by utilizing the secular perturbation theory, 
developed by Mardling (2007). The particular expan- 
sion of the disturbing function in terms of semi-major 
axes ratios places no restriction on the outer eccentricity 
in the equations of motion. Consequently, we can solve 
for the eccentricity ratio of the two planets by explicitly 
equating the precession rates of the two planets, given 
by 
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Fig. 3. — A modified, dissipated Laplace-Lagrange secular solu- 
tion of the 61 Vir system. The initial conditions were identical to 
those, listed in table (1). 

where v = cos(tui — zu 2 ) = ±1. These equations ex- 
plicitly reveal that (ei/e 2 ) is not independent of e 2 , as 
suggested by the eigenvector solutions. Note however, 
that the same expression for the eigenvectors (18)-(19) 
can be derived from these equations by expanding them 
to linear order in e and solving for (ei/e 2 ). The solutions 
for (ei/e 2 ) as a function of /, obtained using equations 
(20) and (21) are also shown in figure (2) for various 
values of e 2 . 

2.4. The General Case: Gauss's Averaging Method 

The above examples are illustrative in nature and are 
applicable when the appropriate assumptions are satis- 
fied. It is also useful, however, to consider a general 
method that will be applicable in all cases, as long as 
the interactions among the planets are secular in nature. 

Rather than expanding the disturbing function in 
terms of a small parameter and applying Lagrange's 
planetary equations, consider N coplanar interacting el- 
liptical wires of mass where the line density is inversely 
proportional to orbital speed and the integrated mass of 
the wire amounts to that of the planet (Gauss 1818). 
The magnitude of the force exerted on line element rjdfj 
by a line element r k df k is simply 



F jk - G^-j°^ k df J df k 
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where r is orbital radius, p is density, / is true anomaly, 
and A = \rj — r k \ is the distance between the line ele- 
ments. The radial and a tangential components of the 
force on line elements j and k are then 
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where 4> = (/fe+^fe - fj~ w j) ls the angle between the line 
elements (Murray & Dermott 1999). Recall that we are 
only interested in the situation where w\ — vj 2 = (0,7r). 
Following Burns (1976), the perturbation equation for 
the precession for longitude of perihelion reads 
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Fig. 4. — Fixed point eccentricity ratio of planet b to planet d as 
a function of the 61 Vir system inclination. The curves were com- 
puted using Gaussian averaging method, with different stationary 
eccentricities of planet d, as labeled. 

with an identical equation for dw^/dt. Note that in this 
formulation, as before, the secular term is linearly pro- 
portional to sin(7), unlike the GR correction. Thus, the 
system inclination can be solved for in the same way as 
above, but without constraints on eccentricity of semi- 
major axes. 

3. APPLICATION: 61VIR 

To date, the number of detected multi-planet systems 
that host small close-in planets remains limited to a 
handful of systems: HD 40307, 55 Cnc, 61 Vir, GJ 581 
and GJ 876. Furthermore, the data for these systems 
are still comparatively sparse, so the error bars on the 
planet's eccentricities are rather large. These issues will 
surely get resolved with time, but at this point we can 
only give a rough assessment, and shall limit our analysis 
to a single case: 61 Vir. 

The planetary system around the nearby sun-like star 
61 Vir was discovered by Vogt et al (2010). The star 
hosts 3 planets, with orbital periods of roughly 4. 2d, 
38d and 124d (see Table 1 for an orbital fit). A sim- 
ple evaluation of the system's dynamical stability yields 
no useful constraints on the inclination of the system. 
However, the minimum mass of the inner-most planet of 
m = 5.1 ± 0.6M e corresponds to that of a super-Earth, 
making it an ideal candidate for our method. 

The characteristic isolated circularization timescale of 
planet b is roughly r c ~ Q x 10 6 years. A damped, 
modified LL solution (shown in Figure 3) reveals that 
depending on starting conditions, up to 10 r is required 
for the system to arrive to the fixed point. Thus, as al- 
ready pointed out by Vogt et al (2010), given the star's 
multi-billion year age, we expect the system to be sta- 
tionary if Qf, < 10 3 . For the illustrative purposes of this 
paper, we assume that planet b's tidal quality factor is 
similar to that of rocky bodies i.e. Qb = 100. 

Initially, we proceed as described in section 2.1 and 
compute the surviving LL eigenvector that physically 
corresponds to a state where all orbits are apsidally 
aligned. Given the moderate eccentricity (e > 0.1) of 
the outer two planets, however, the LL solution does not 
give a quantitatively acceptable answer. Consequently, 
we recompute the eccentricity ratios using the Gaussian 
averaging method, as described in section 2.4, utilizing 
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Fig. 5. — Fixed point eccentricity ratio of planet c to planet d as 
a function of the 61 Vir system inclination. The curves were com- 
puted using Gaussian averaging method, with different stationary 
eccentricities of planet d, as labeled. 



TABLE 1 
Orbital Fit of the 61 Vir System 



Planet Mass (019) P (days) 



(deg) 



5.28 
19.1 
23.4 



4.3 

38 
123 



0.147 
0.155 
0.34 



104 

331 
314 



the LL solution as an initial guess in the root-finding al- 
gorithm. The resulting curves are plotted in figures (4) 
and (5). It is noteworthy that although the Gaussian 
and LL solutions are qualitatively similar, higher-order 
secular terms clearly make a noticeable contribution to 
the fixed-point solution. 

Although the error bars on the orbital elements are 
still large, it is noteworthy that the observed system is 
consistent with a fixed point configuration. Thus, further 
observation of the system is warranted, given that if the 
system is found to be in a stationary state, it would yield 
not only the true masses, but also a constraint on the 
tidal quality factor of the inner- most planet. 



4. COMMENTS ON MASSIVE PLANETS 

The domain of applicability of the method described in 
this paper does not extend to "large" planets (recall that 
we require A^ idal <C 1 in order to solve for sin(/)). How- 
ever, for massive, close-in planets, the sin(7) degeneracy 
can be resolved from spectral characterization of the host 
star alone (Snellen et al 2010). In such a case, the orbital 
precession rate yields information on the radius and the 
interior structure of the planet. 

If only a single planet is present in the system, then 
the method described by Ragozzine & Wolf (2009) can 
be employed. Namely, if the planet is sufficiently close to 
its host star, the orbital precession rate may be as high 
as a few degrees/year. In this case, direct observation 
of the orbital precession can be related to the sum of 
equations (7) -(9). As already discussed above, however, 
the first term in equation (9) dominates all other terms 
for large, massive planets. Consequently, k2(R) 5 can be 
inferred. 

In order to accurately measure orbital precession, es- 
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Fig. 6. — Parameter space over which the method described here 
is applicable. The purple region is characterized by planets with 
circularization timescales less that 1 Gyr and non-relativistic ef- 
fects contributing to less than 10% of the non-secular precession. 
The blue dots correspond to currently known low-mass RV planets. 
Contours of RV signal semi-amplitudes are also shown. 

pecially within the context of RV observations, signifi- 
cantly non-zero orbital eccentricity is needed. This poses 
a problem, since the eccentricities of single close-in plan- 
ets are usually damped out on the timescale of ~ 1 Gyr. 
As a result, in practice, the method of Ragozzine & Wolf 
(2009) is much better suited for transiting planets, where 
ultra-precise photometry, such as that characteristic of 
the Kepler mission, can be used to pinpoint even a low 
(~ 10~ 3 ) eccentricity. 

If there are two or more planets in the system, the 
situation is considerably more advantageous, since a fi- 
nite eccentricity of the inner planet can be maintained 
over the age of the star by a perturbing companion. In 
this case, under the assumption of co-planar planets, the 
characterization of the fixed point through equations (20) 
and (21), where B\\ is dominated by the planetary tidal 
term, yields /^(-R) 5 - In essence, the calculation is analo- 
gous to that of Batygin, Bodenheimer & Laughlin (2009) 
for the Hat-P-13 system, with the exception that the ra- 
dius is also unknown. Unfortunately, for a given mass, 
k 2 (R) 5 is not a single- valued function of R so the values 
of hi and R cannot be disentangled by modeling of the 
planetary interior. 

Determination of k 2 R 5 of interest because the num- 
ber of RV systems where the calculation is applicable is 
bound to greatly exceed the number of transiting systems 
for which fc 2 can be measured directly, and a substantial 
distribution can be formed. The results of the Kepler 
mission will provide a statistical distribution for plan- 
etary radii. However, because the majority of stars in 
the Kepler field of view are faint, RV follow-up of most 
systems will be difficult. This poses a challenge for de- 
termination of &2 by the method proposed by Batygin, 
Bodenheimer & Laughlin (2009). Consequently, there is 
considerable value in deriving a statistical distribution 
for Love numbers from these observations. 

5. CONCLUSION 

In this paper, we present a method for determination of 
the true masses of RV planetary systems with a close-in 
planet. The analysis in question has important implica- 
tions. First and foremost, it opens up a direct avenue 
towards an identification of the true lowest-mass exo- 



planets detected. This provides a direct constraint on 
the discussion of the habitability of RV planets. The 
second implication is more indirect. In a recent study, 
Ho & Turner (2010) showed that there is significantly 
more uncertainty in sin(/) than previously assumed. In 
particular, the assumption that sin(7) has been drawn 
from a flat distribution is incorrect. Instead the distribu- 
tion from which sin(J) is drawn is sensitively dependent 
on the true masses. Without additional information, it 
appears that there is significant adversity in estimating 
true masses of exoplanets from observations alone. As 
a result, resolution of sin(/) using an independent tech- 
nique holds value not only in just yielding the true masses 
of a particular system, but also in implicitly constraining 
the relevant true-mass distribution from which sin(J) can 
then be drawn for the entire sample. 

It is certainly worthwhile to consider the observability 
of the systems to which our method is applicable. Recall 
that our method relies on three assumptions. First, tidal 
dissipation of orbital energy by the inner-most planet 
of a system must be efficient enough for the system to 
become tidally relaxed on a time-scale, less than a few 
Gyr i.e. the age of the star. Second, we require rough 
coplanaritjj^] of the system to ensure that fixed-point ec- 
centricities are unaffected by the precession of the as- 
cending node (Mardling 2010). Finally, to separate the 
dependence on sin(/) in the equations of motion, we re- 
quire that the additional precession of the perihelion of 
the inner-most planet arises primarily from GR. Upon 
satisfaction of the above criteria, sin(7) can be solved for 
in an explicit, direct way. 

To demonstrate the extent of parameter space over 
which our method is applicable, we delineated the re- 
gion where non-GR contributions account for less than 
10% of the additional precession of the inner-most planet, 
and circularization timescale is less than 1 Gyr. Figure 
(6) shows this range, along with the current aggregate of 
low-mass RV planets. Given the uncertainty in tidal Q 
as a function of planetary mass, r c = 2 Gyr and r c = 0.5 
Gyr curves are also presented. Additionally, contours 
of corresponding semi-amplitudes of RV signal (K) are 
also displayed. Although the parameter space covered is 
considerable, it is clear that approximately 3-day period 
hot Neptunes make the best candidates for our method 
because of the optimum interplay between K (making 
the planets most readily observable), and r c . Finally, 
we discuss the possibility of obtaining information about 
the radius and interior structure of massive hot Jupiters 
in multiple systems, where the sin(7) degeneracy can be 
resolved with observations alone. Consequently, we en- 
courage continued RV observation and more importantly, 
follow-up of qualifying multi-planet systems, with the 
goal to pinpoint the orbital state to a high precision, 
thus deriving true masses and constraining the interior 
structure of low and high-mass RV exoplanets, respec- 
tively. 

We would like to conclude by presenting a list of 
possibilities for determination of physical properties of 
planets from observations of orbital parameters. The 
compiled flow-chart is presented as Figure 7. Let us 

2 "Rough" coplanarity implies that the terms in the disturbing 
function that have the mutual inclination as a multiplier are small 
in comparison with terms of the same order that contain only the 
eccentricities. 
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summarize: if a newly discovered system harbors only 
a single tansiting hot Jupiter, the interior structure 
can be derived from monitoring of orbital precession 
(Ragozzine & Wolf 2009). Alternatively, although 
observationally challenging, the rotational and tidal 
bulges can be deduced directly from the shape of the 
light-curve (Carter & Winn 2010, Leconte et al 2011). If 
two planets are present and reside at a fixed point, the 
situation becomes more advantageous. If tidal precession 
plays an important role, and the inner planet transits, 
the Love number can be derived from a single snap-shot 
observation of the orbital state (Batygin Bodenheimer & 
Laughlin 2009). If the inner planet does not transit, its 
exact mass can be derived spectroscopically (Snellen et 
al 2010) and fc2-R 5 can be computed. On the other hand, 
if GR overwhelms tidal precession, sin(J) degeneracy 
of the system can be resolved. If the system is tidally 
relaxed but is not co-planar, orbital evolution will follow 
a limit cycle rather than a fixed point (Mardling 2010). 
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Fig. 7. — A flow-chart that depicts various possibilities for de- 
termination of physical properties of planets from observations of 
their orbital configurations. 



